IMPLANTABLE DEVICES AND METHODS USING FREQUENCY-DOMAIN 
ANALYSIS OF THORACIC SIGNAL 



TECHNICAL FIELD 

5 This document relates generally to medical systems, devices, and methods, 

and particularly, but not byway of limitation, to implantable devices and methods 
using frequency-domain analysis of a thoracic signal. 

BACKGROUND 

10 Implantable cardiac rhythm management devices typically monitor and 

process heart signals to provide the heart with needed therapy. One example of such 
therapy includes delivering a pacing pulse that triggers a resulting heart contraction. 
Another example includes delivering a defibrillation shock that interrupts an 
abnormal heart rhythm. Some cardiac rhythm management devices also monitor the 

15 patient's breathing (also referred to as the patient's respiration). Such devices 
typically use breathing to control the delivery of pacing pulses to increase the 
patient's heart rate when the patient breathes more rapidly, and to slow the patient's 
heart rate when the patient breathes more slowly. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

In the drawings, which are offered by way of example, and not by way of 
limitation, and which are not necessarily drawn to scale, like numerals describe 
substantially similar components throughout the several views. Like numerals 
having different letter suffixes represent different instances of substantially similar 
25 components. 

FIG. 1 is a schematic diagram illustrating generally one example of portions 
of an implantable unit, such as implantable cardiac rhythm management device. 

FIG. 2 is a schematic diagram illustrating generally further details of portions 
of one example of the implantable device. 
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FIGS. 3A, 3B, and 3C jointly provide a generalized conceptual illustration of 
normal breathing. 

FIGS. 4A, 4B, and 4C jointly provide a generalized conceptual illustration of 
periodic breathing. 

5 FIGS; 5 A, 5B, and 5C jointly provide a generalized conceptual illustration of 

Cheyne-Stokes breathing. 

FIG. 6 is a flow chart illustrating generally one technique for using the 
implantable device for classifying a subject's respiration state. 

FIG. 7 is a flow chart illustrating generally one technique for analyzing the 
10 respiration spectrum. 

FIG. 8 is another generalized conceptual illustration of normal breathing. 
FIG. 9 is a generalized conceptual illustration of an abnormal respiration 
state, which exhibits shortened inhalation and prolonged exhalation. 

FIG. 10 is a flow chart illustrating generally one technique of performing 
15 frequency domain spectral analysis on the respiration spectral signal, such as to 

distinguish the abnormal respiration state of FIG. 9 from the normal respiration state 
of FIG. 8. 

FIG. 11 is a conceptualized graph of a transthoracic impedance spectral 
signal amplitude vs. frequency for a resting subject. 
20 FIG. 12 is a conceptualized graph of a transthoracic impedance spectral 

signal amplitude vs. frequency for an exercising subject. 

FIG. 13 A is a generalized conceptual illustration of a transthoracic 
impedance signal in the time and frequency domains for a subject without fluid 
accumulation in or around the lungs. 
25 FIG. 13B is a generalized conceptual illustrations of a transthoracic 

impedance signal in the time and frequency domains for a subject with fluid 
accumulation in or around the lungs. 

FIG. 14 is a schematic diagram, similar to FIG. 2, but illustrating a thoracic 
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monitor circuit that obtains a thoracic signal from an acceleration detector circuit, 

DETAILED DESCRIPTION 
In the following detailed description, reference is made to the accompanying 
5 drawings, which form a part hereof, and in which is shown by way of illustration 
specific embodiments in which the invention may be practiced. These 
embodiments, which are also referred as "examples," are described in sufficient 
detail to enable those skilled in the art to practice the invention, and it is to be 
understood that the embodiments may be combined, or that other embodiments may 

10 be utilized and that structural, logical and electrical changes may be made without 
departing from the scope of the present invention. The following detailed 
description is, therefore, not to be taken in a limiting sense, and the scope of the 
present invention is defined by the appended claims and their equivalents. 

In this document, the terms "a" or "an" are used, as is common in patent 

15 documents, to include one or more than one. Furthermore, all publications, patents, 
and patent documents referred to in this document are incorporated by reference 
herein in their entirety, as though individually incorporated by reference. In the 
event of inconsistent usages between this documents and those documents so 
incorporated by reference, the usage in the incorporated reference(s) should be 

20 considered supplementary to that of this document; for irreconcilable 
inconsistencies, the usage in this document controls. 

This document discusses, among other things, systems, devices, and methods 
that will be described in applications involving implantable medical devices 
including, but not limited to, implantable cardiac rhythm management devices such 

25 as pacemakers, cardioverter/defibrillators, pacer/defibrillators, biventricular or other 
multi-site ^synchronization or coordination devices, and drug delivery systems. 
However, these systems, devices, and methods may be employed in unimplanted 
devices, including, but not limited to, external pacemakers, cardioverter/ 
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defibrillators, pacer/defibrillators, biventricular or other multi-site ^synchronization 
or coordination devices, monitors, programmers and recorders, whether such 
devices are used for providing a diagnostic, a therapy, or both a diagnostic and a 
therapy. 

5 The present inventors have recognized that certain sicknesses also affect 

pulmonary function or both pulmonary and cardiac function, and that there exists an 
unmet need for providing improved systems, devices, and methods for using 
information obtained from the patient's pulmonary function or both pulmonary and 
cardiac function to provide better diagnosis and/or treatment of certain patients 

1 0 and/or for determining therapy efficacy. 

FIG. 1 is a schematic diagram illustrating generally, by way of example, but 
not by way of limitation, portions of an implantable unit, such as implantable 
cardiac rhythm management device 100. The example of FIG. 1 also illustrates 
portions of one environment in which the device 100 is capable of being used. In 

1 5 the example illustrated in FIG. 1, the implantable unit 100 is operatively coupled to 
the subject's heart 102 and to a portion of the subject's thorax using implantable 
electrodes 104. The electrodes 104 are located with the device 100 and/or are 
located in association with the heart, such as by being disposed on the distal end of 
an intravascular leadwire 106 that is coupled to the device 100. The example of 

20 FIG. 1 also illustrates an external remote interface unit 108. In this example, the 
remote interface unit 108 is located external to the subject. The remote interface 
unit 108 is wirelessly or otherwise communicatively coupled to a telemetry circuit 
110 included in the implantable device 100. 

In the example of FIG. 1, the device 100 includes a thoracic monitor circuit 

25 112. As discussed below, one example of the thoracic monitor circuit 112 includes 
an impedance detector circuit. Another example of the thoracic monitor circuit 112 
includes an acceleration detector circuit. In the example in which the thoracic 
monitor circuit 112 includes an impedance detector circuit, the impedance detector 
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circuit is configured to use the implanted electrodes 104 to detect a transthoracic 
impedance associated with a portion of the subject's thorax. One suitable electrode 
configuration for detecting transthoracic impedance is described in commonly- 
assigned Hartley et al. U.S. Patent Number 6,076,015 entitled RATE ADAPTIVE 
5 CARDIAC RHYTHM MANAGEMENT DEVICE USING TRANSTHORACIC 

IMPEDANCE, which is incorporated herein by reference in its entirety, including its 
description of an electrode configuration and circuitry to detect transthoracic 
impedance. Another suitable electrode configuration for detecting transthoracic 
impedance is described in commonly-assigned Hauck et al. U.S. Patent No. 

1 0 5,284, 136 entitled DUAL INDIFFERENT ELECTRODE PACEMAKER, which is 
incorporated herein by reference in its entirety, including its description of an 
electrode configuration and circuitry to detect transthoracic impedance. The 
electrode configuration need not use a header electrode 104A. In different 
examples, the electrode configuration includes two electrodes 104, three electrodes 

15 104, and four electrodes 104. 

The thoracic monitor circuit 112 includes an output 114 to provide a thoracic 
signal. The thoracic signal provided at the output 1 14 includes information about 
the subject's thoracic impedance and/or acceleration. Such thoracic information 
provides information about the subject. For example, one component of such 

20 thoracic information provides information about the subject's heart contractions, 
including information about the rate at which such atrial and/or ventricular heart 
contractions occur. Another component of such thoracic information provides 
information about the subject's breathing or respiration. 

In the example of FIG. 1, the device 100 also includes a digital signal 

25 processor (DSP) circuit 116. The DSP circuit 116 includes an input 118 that is 

coupled to the output 114 of the thoracic monitor circuit 112 to receive a digitized 
thoracic signal provided by the thoracic monitor circuit 112. As discussed below, the 
DSP circuit 116 performs frequency domain signal processing on the thoracic signal. 
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One example of such frequency domain signal processing includes separating 
cardiac contraction information from lower frequency components of the thoracic 
signal. Another example of such frequency domain signal processing includes 
extracting respiration information from the thoracic signal In a further example, the 
5 DSP circuit 116 uses spectral analysis to classify the subject's pulmonary 

physiological state. The DSP circuit 116 includes an output 120 that provides 
information about the subject, for example, information about the pulmonary 
physiological classification. In a further example, the telemetry circuit 110 includes 
an input 122 that is coupled to the output 120 of the DSP circuit 116 to receive an 

1 0 indication of the pulmonary physiological classification or the like. In this example, 
the telemetry circuit 110 is capable of communicating at least one indication of the 
pulmonary physiological classification to the external remote interface unit 108, for 
external storage and/or display, either of a single such indication instance or of a 
trend of such indication instances. 

15 FIG. 2 is a schematic diagram illustrating generally, by way of example, but 

not by way of limitation, further details of portions of one example of the 
implantable device 100. In the example of FIG. 2, the thoracic monitor circuit 112 
includes an impedance detecting circuit 201, comprising a test stimulus circuit 200 
to provide an excitation or stimulus to a portion of the subject's thorax for detecting 

20 transthoracic impedance. The impedance detecting circuit 201 also includes a 
response sensing circuit 202 to sense a responsive indication correlative to 
transthoracic impedance. In one example, the test stimulus circuit 200 includes a 
test current carrier signal generator, and the response sensing circuit 202 includes a 
voltage preamplifier circuit, such as described in detail in Hartley et al. U.S. Patent 

25 Number 6,076,015, which is incorporated by reference herein in its entirety, 
including its description of an exciter circuit and a preamplifier circuit. 

In the example illustrated in FIG. 2, at least one output 204 of the response 
sensing circuit 202 is coupled to at least one input 206 of an analog signal 
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processing circuit 208. In one example, the analog signal processing circuit 208 
includes, among other things, a demodulator circuit and a bandpass filter circuit, 
such as described in detail in Hartley et al. U.S. Patent Number 6,076,015, which is 
incorporated by reference herein in its entirety, including its description of 
5 demodulator and bandpass filter circuits. In the example of FIG. 2, an output 210 of 
the analog signal processing circuit 208 is coupled to an input 212 of an analog-to- 
digital (A/D) converter circuit 214, such as the 8-bit A/D converter described in 
Hartley et al. U.S. Patent Number 6,076,015, which is incorporated by reference 
herein in its entirety, including its description of an A/D converter circuit. At least 

10 one output 216 of the A/D converter circuit 214 is coupled to the output 114 of the 
thoracic monitor circuit 112 to provide a digitized thoracic signal coupled by a bus 
connection or otherwise to the input 118 of the DSP circuit 116. 

In the example of FIG. 2, the DSP circuit 116 includes a time-domain-to- 
frequency-domain conversion circuit, such as fast Fourier transform (FFT) module 

15 218. The FFT module 218 may be implemented either as a dedicated hardware 

circuit or as a sequence of instructions executed on a microprocessor circuit or the 
like. The FFT module 218 includes an input 220 that is coupled to the input 118 of 
the DSP circuit 116 to receive the digitized time-domain thoracic signal provided by 
the thoracic monitor circuit 112. In this example, the FFT module 218 converts the 

20 digitized time-domain thoracic signal received from the thoracic monitor circuit 112 
into a digitized frequency-domain thoracic signal provided at an output 222 of the 
FFT module 218. 

In the example of FIG. 2, the DSP circuit 116 includes a frequency-domain 
filter 224. The filter 224 may be implemented either as a dedicated hardware circuit 
25 or as a sequence of instructions executed on a microprocessor circuit or the like. An 
input 226 of the frequency-domain filter 224 is coupled to the output 222 of the FFT 
module 218 to receive the frequency-domain thoracic signal. In one example, the 
frequency-domain filter 224 extracts a frequency-domain respiration signal from the 
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frequency-domain thoracic signal (which may also include confounding signal 
content, such as a higher frequency heart contraction component, and such as a 
lower frequency component due to fluid shifts toward or away from the thorax). 
The filter 224 provides the resulting frequency-domain respiration signal at an 
5 output 228. 

In the example of FIG. 2, the frequency domain filter 224 is implemented as 
an adaptive frequency-domain filter 224 in which a filter cutoff frequency varies 
according to a detected heart rate. In this example, the device 100 includes at least 
one sense amplifier or other such depolarization detector circuit 230. The 

10 depolarization detector circuit 230 includes inputs 232 that are coupled to the 

implanted electrodes to detect heart depolarizations indicative of heart contractions. 
The depolarization detector circuit 230 includes one or more outputs 234 providing 
interrupts triggered by such detected heart depolarizations. A heart rate interval 
timer circuit 236 includes one or more corresponding inputs 238 coupled to the 

1 5 respective outputs 234 of the depolarization detector circuit 230 to receive the 
interrupts. The heart rate interval timer circuit 236 computes the time interval 
between interrupts indicating heart depolarizations. An output 239 of the heart rate 
interval timer circuit 236 provides information correlative to the heart rate frequency 
between successive heart contractions to an input 240 of the frequency-domain 

20 adaptive filter 224. In response to such heart rate information, the frequency- 
domain adaptive filter 224 adjusts a cutoff frequency in the frequency domain. In 
one example, this heart-rate dependent variable cutoff frequency is used for 
separating, in the frequency domain, heart contraction information from respiration 
information, such as to obtain a frequency-domain signal correlative to the subject's 

25 heart contractions and respiration, respectively. In this example, heart contraction 
information is provided at output 229 of the filter 224. Respiration information is 
provided at the output 228 of the filter 224. 
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In the example of FIG. 2, the DSP circuit 116 includes a frequency-domain 
spectrum analyzer 242. The spectrum analyzer 242 may be implemented either as a 
dedicated hardware circuit or as a sequence of instructions executed on a 
microprocessor circuit or the like. The spectrum analyzer 242 includes an input 244 
5 that is coupled to receive the frequency domain respiration signal from the output 
228 of the filter 224. In one example, the spectrum analyzer 242 classifies the 
subject's respiration using frequency domain spectrum analysis of the frequency 
domain respiration signal. An output 246 of the spectrum analyzer 242 provides a 
resulting digitally encoded pulmonary physiological state classification signal (e.g., 

10 normal breathing, periodic breathing, Cheyne-Stokes breathing, obstructed 

respiration, restrictive respiration, pulmonary fluid accumulation, etc.) In one 
example, this pulmonary physiological classification is communicated to the 
telemetry circuit 110, for further communication from the subject to the external 
remote interface unit, such as for assisting a caregiver in diagnosing at least one 

15 aspect of the subject's health status. In another example, the pulmonary 

physiological classification is communicated to an input 248 of a controller 250. 
The controller 250 includes one or more outputs 252 that are coupled to 
corresponding inputs 254 of a therapy circuit 256, such as for controlling its delivery 
of pacing, defibrillation, ^synchronization or other cardiac rhythm management 

20 therapy provided to the subject using implanted electrodes. In one example, the 

controller 250 adjusts the cardiac rhythm management therapy using, among other 
things, the pulmonary physiological classification information. 

In one further variation of the above system 100, the DSP circuit 116 also 
includes a frequency domain heart rate variability (HRV) spectrum analyzer circuit 

25 258 for processing heart contraction information received at an input 257 from the 
output 229 of the frequency domain adaptive filter 224. This embodiment is 
discussed further below. 
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FIGS. 3A, 3B, and 3C jointly provide a generalized conceptual illustration of 
normal breathing. FIG. 3 A illustrates a conceptual example of a time domain 
normal respiration signal, illustrating respiration amplitude vs. time. FIGS. 3B and 
3C provide corresponding spectral illustrations (amplitude vs. frequency) of this 
5 normal respiration signal after being converted into the frequency domain. In this 

example, normal respiration is characterized by a respiration rate of between about 5 
breaths per minute and about 20 breaths per minute, when the subject is at rest. 
However, when the subject is exercising, the respiration rate will typically be 
between about 15 breaths per minute and about 60 breaths per minute. FIG. 3B 

10 illustrates an idealized frequency spectrum, for the resting subject, having the energy 
concentrated at about 0.17 Hz (i.e., about 10 breaths per minute), and assumes pure 
sinusoidal respiration at a fixed rate. FIG. 3C illustrates the frequency spectrum of a 
more likely encountered respiration signal that is not a pure sinusoid and includes 
some variation in the respiratory rate. 

15 FIGS. 4 A, 4B, and 4C jointly provide a generalized conceptual illustration of 

periodic breathing. FIG. 4A illustrates a conceptual example of a time-domain 
periodic respiration signal, illustrating respiration amplitude vs. time. FIGS. 4B and 
4C provide corresponding spectral illustrations (amplitude vs. frequency) of this 
periodic respiration signal after being converted into the frequency domain. In this 

20 example, periodic respiration is characterized by a modulation in breathing 

amplitude. The modulation typically exhibits a period between about 0.5 minutes 
and about 2 minutes. FIG. 4B illustrates the idealized resulting frequency spectrum, 
for the resting subject, having most energy concentrated at about 0,17 Hz (i.e., about 
10 breaths per minute), but also having considerable energy concentrated at about 

25 0.167 Hz - 0.017 Hz = 0.150 Hz and at about 0.167 Hz + 0.017 Hz - 0.184 Hz (i.e., 
at the breathing frequency of about 1 cycle per minute, corresponding to the 
modulation frequency). 



Attorney Docket No. 279.580US1 



} 0 



FIGS. 5 A, 5B, and 5C jointly provide a generalized conceptual illustration of 
Cheyne-Stokes breathing. FIG. 5A illustrates a conceptual example of a time- 
domain Cheyne-Stokes respiration signal, illustrating respiration amplitude vs. time. 
FIGS. 5B and 5C also illustrate a corresponding spectral illustration (amplitude vs. 
5 frequency) of this Cheyne-Stokes respiration signal after being converted into the 

frequency domain. In this example, Cheyne-Stokes respiration is characterized by a 
modulation in breathing amplitude, similar to that illustrated in FIG. 4 for periodic 
breathing. The modulation typically exhibits a period of between about 0.5 minutes 
and about 2 minutes. However, unlike periodic breathing, the modulation that 

10 occurs in Cheyne-Stokes breathing results in the respiration amplitude going to zero 
during the "valleys" of the modulated respiration signal, as illustrated in FIG. 5A. 
FIGS. 5B and 5C illustrate the resulting frequency spectrum for the resting subject. 
Like the case of periodic breathing, most of the spectral energy illustrated in 
Cheyne-Stokes breathing of FIGS. 5A and 5B is concentrated at about 0.17 Hz (i.e., 

15 about 10 breaths per minute). Moreover, the Cheyne-Stokes breathing of FIGS. 5B 
and 5C includes even less energy concentrated at 0.167 Hz than for the case of 
periodic breathing. 

The idealized spectrums shown in FIG. 4B and FIG. SB represent a 
respiration signal that is purely sinusoidal and at a constant rate. FIG. 4B and FIG. 

20 5B also represent a modulation signal (i.e., the periodic breathing) that is purely 

sinusoidal, where the modulation of the respiratory signal by periodic breathing can 
be represented by a purely multiplicative process. FIG. 4C and FIG. 5C illustrate 
the frequency spectrum of a more likely encountered periodic respiration and 
Cheyne-Stokes respiration, respectively. FIG. 4C and FIG. 5C represent modulation 

25 signals that are not pure sinusoids and also assumes some variation in respiratory 

rate. In addition, energy is shown in FIG. 4C and FIG. 5C at the periodic breathing 
frequency (i.e., about 1 cycle per minute) due to an additive modulation of the 
respiration signal by the periodic breathing signal. Finally, FIG. 4C and FIG. 5C 
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illustrate a merging of the numerous frequency components around 0.167 Hz due to 
finite resolution of the spectral analysis apparatus. 

FIG. 6 is a flow chart illustrating generally, by way of example, but not by 
way of limitation, one technique for using the device 100 for classifying a subject's 
5 respiration state. At 600, a thoracic signal is acquired by the thoracic monitor circuit 
112, such as by using the techniques that are described in detail in Hartley et al. U.S. 
Patent Number 6,076,015, which is incorporated by reference herein in its entirety, 
including its description of techniques for acquiring a transthoracic impedance 
signal. At 602, the thoracic signal is digitized, such as by the A/D converter circuit 

10 214. At 604, the digitized time-domain thoracic signal is converted into a 

frequency-domain thoracic spectral signal, such as by using the FFT module 218. 
At 608, a frequency domain filtering is performed to extract the respiration spectral 
components of the thoracic spectrum. At 610, this resulting respiration spectrum is 
analyzed in the frequency-domain, such as to obtain information about the subject's 

15 respiration state. At 612, the respiration state is classified using information 
obtained from the frequency-domain analysis of the respiration spectrum. 

One consideration in performing the time-domain-to-frequency-domain 
conversion of 604 is in selecting the appropriate time window of the time domain 
thoracic signal for converting into a frequency domain thoracic spectral signal (e.g., 

20 by performing an FFT) for subsequent analysis of respiration spectral information. 
In one example, a short term analysis of respiration is desired (e.g., a time window 
that includes several breaths, and up to about 1 minute, or several minutes). This 
example permits determining the pulmonary physiological state of the subject. For 
example, it allows discrimination between normal breathing, periodic breathing, and 

25 Cheyne-Stokes breathing, etc. In another example, however, a longer term analysis 
of respiration is desired. In one such example, the time window may include days or 
even months. This longer term analysis of respiration may be used to determine the 
pulmonary physiological state of the subject, such as described above, over the 
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longer term period. Alternatively, such longer term analysis may be used to provide 
a more generalized indication of the patient's health status. In one such example of 
a generalized patient health status indicator, a wide respiration spectral dispersion 
with a predominance of low frequency energy is used as an indicator of a relatively 
5 healthy state, while a narrow respiration spectral dispersion with relatively high 

frequency energy is used as an indicator of a relatively unhealthy state. In another 
example, substantially continuous high rate breathing, particularly with a low tidal 
volume (amplitude), will provide a sensitive indication of pulmonary fluid 
accumulation, also referred to as pulmonary edema. 

10 FIG. 7 is a flow chart illustrating generally, by way of example, but not by 

way of limitation, one technique for analyzing the respiration spectrum, such as 
illustrated at 610 of FIG. 6. In FIG. 7, at 700, a peak amplitude of the frequency- 
domain respiration spectral signal is obtained. In one example, the frequency 
corresponding to this peak amplitude is also obtained. At 702, modulation band test 

15 thresholds 77 and T2 are obtained using the spectrum analyzer 242. In one example, 
the test thresholds Tl and T2 are spectral amplitudes that are normalized (e.g., a 
specified percentage) to the previously-obtained peak amplitude of the respiration 
spectrum. In another example, the test thresholds Tl and T2 are energy thresholds 
. that are normalized (e.g., a specified percentage) to a spectral energy that is obtained 

20 by actually or approximately integrating a spectral energy density over a spectral 
interval [fl 9 /2] about the frequency (f 0 ) of the peak amplitude. 

At 704, the modulation band (e.g., breathing at frequencies having 
corresponding periods between about 0.5 minutes and about 2 minutes, as discussed 
above) spectral signal is examined. For the case in which the test thresholds Tl and 

25 T2 include spectral amplitudes, then the examination of the modulation band at 704 
includes determining a peak amplitude in the modulation band. For the case in 
which the test thresholds Tl and T2 include energy thresholds, then the examination 
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of the modulation band at 704 includes integrating the spectral energy across the 
frequencies in the spectral interval \f3,f4] of the modulation band. 

At 706, the modulation band signal is compared to the first threshold 77. If 
the modulation band peak amplitude is less than the amplitude threshold 77 (or, 
5 alternatively, the integrated spectral energy in the modulation band is less than the 
energy threshold 77), then at 708, the respiration state is classified as normal 
respiration. 

Otherwise, at 710, the modulation band signal is compared to the second 
threshold T2 (which has a greater magnitude than the threshold 77). If the 
1 0 modulation band peak amplitude exceeds the amplitude threshold T2 (or, 

alternatively, the integrated spectral energy in the modulation band exceeds the 
energy threshold 77), then, at 712, the respiration state is classified as Cheyne- 
Stokes respiration. Otherwise, at 714, the respiration state is classified as periodic 
breathing. 

1 5 FIG. 8 is a generalized conceptual illustration of normal breathing, similar to 

FIG. 3, but illustrating a higher time-domain resolution. FIG. 9 is a generalized 
conceptual illustration of another abnormal respiration state (with a time-domain 
resolution comparable to FIG. 8). The abnormal respiration state of FIG. 9 exhibits 
decreased inhalation time (a steeper slope on the rising "inhalation" portion of the 

20 transthoracic impedance signal), and an increased exhalation time (a more gradual 
descent on the falling "exhalation" portion of the transthoracic impedance signal). 
Acutely, this may result from an asthma episode. Chronically, this may result from 
pulmonary obstructive disease (such as asthma, emphysema, and/or bronchitis). The 
time-domain abnormal respiration signal illustrated in FIG. 9 results in harmonic 

25 content occurring at a higher "harmonic" sideband frequency than the peak 

amplitude frequency f 0 of about 0.17 Hz. By contrast, a restrictive respiration 
disease (such as pectus excavatum, myasthenia gravis, diffuse idiopathic interstitial 
fibrosis, and space occupying lesions such as tumors and/or effusions) is 
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characterized by a reduced tidal volume, which results in a reduced spectral 
amplitude than for normal breathing, and faster breathing, which shifts the spectral 
content of the respiration signal to a higher frequency than for normal breathing. 

FIG. 10 is a flow chart illustrating generally, by way of example, but not by 
5 way of limitation, one technique of performing frequency domain spectral analysis 
on the respiration spectral signal, such as to distinguish the abnormal respiration 
state of FIG. 9 from the normal respiration state of FIG. 8. In the example 
illustrated in FIG. 10, at 1000, the peak amplitude of the respiration spectral signal is 
obtained. In one example, the corresponding frequency of such peak amplitude is 

10 also obtained. At 1002, a harmonic band energy threshold T3 is obtained using the 
spectrum analyzer 242. In one example, the test threshold T3 is an energy threshold 
that is normalized (e.g., a specified percentage) to a spectral energy obtained by 
actually or approximately integrating a spectral energy density over a spectral 
interval \fl 9 f2] about the frequency corresponding to the previously-obtained peak 

1 5 amplitude. At 1004, the harmonic band spectral energy is examined. In one 

example, the harmonic band spectral energy is obtained by integrating a respiration 
spectral energy density from a lower frequency limit fS (e.g., approximately equal to 
f2) to a higher frequency limit f6 (e.g., that is approximately equal to the variable 
lowpass cutoff frequency of the frequency-domain adaptive filter 224). At 1006, the 

20 harmonic band spectral energy is compared to the energy threshold T3. If the 
harmonic band spectral energy exceeds T3, then at 1008 the respiration state is 
classified as being abnormal (e.g., indicative of the shortened inhalation, prolonged 
exhalation illustrated in FIG. 9). Otherwise, in one example, at 1010, process flow 
continues at 702 of the flow chart illustrated in FIG. 7, for further classifying the 

25 respiration state, as discussed above. 

FIG. 11 and FIG. 12 are conceptualized graphs of a thoracic spectral signal 
amplitude vs. frequency for a resting subject and an exercising subject, respectively. 
As seen by comparing FIG. 11 with FIG. 12, during exercise, the respiration spectral 
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content moves to a higher frequency that is comparable to the frequency of the heart 
contraction spectral content when the subject is resting. In one example, frequency 
domain adaptive filter 224 of FIG. 2 includes a variable cutoff frequency f c that is 
adjusted as a function of heart rate (HR). The heart rate is obtained from sensed 
5 cardiac depolarizations or, for paced depolarizations, from a therapy circuit control 
signal at 252 that triggers issuance of a contraction-evoking pacing pulse. In one 
example, the variable cutoff frequency is determined by f c = HR/60 - A, where A is a 
fixed frequency offset amount and heart rate is expressed in beats per minute. Heart 
contractions and other spectral content above f c is, in one example, discarded as 

10 noise to the desired respiration spectral content of the thoracic spectral signal. 

However, in an alternative example, the heart contraction spectral content is 
not discarded. Instead, such heart contraction spectral content is provided at the 
output 229 of the frequency domain adaptive filter to the input 257 of the frequency 
domain heart rate spectrum analyzer module 258, as illustrated in FIG. 2. The heart 

15 rate spectrum analyzer module 258 performs heart rate variability (HRV) diagnostic 
testing to provide at least one indication of the subject's health status. Examples 
using spectral analysis of heart beat interval spectral signal information to provide 
HRV diagnostic information are described in commonly-assigned Spinelli et al. U.S. 
Patent No. 5,466,245, Heemels et al. U.S. Patent No. 5,603,33 1, Carlson et al. U.S. 

20 Patent No. 6,026,320, and Carlson et al. U.S. Patent No. 6,301,499, each of which is 
incorporated herein by reference in its entirety, including their disclosed techniques 
of using HRV information, such as to provide an indication of a subject's health 
status. 

FIGS. 13 A - 13B are a generalized conceptual illustrations of a thoracic 
25 impedance signal in the time and frequency domains. FIG. 13A illustrates a time 

domain thoracic impedance signal 1300 and a frequency domain thoracic impedance 
signal 1302 in a subject without pulmonary fluid accumulation in (e.g., pulmonary 
edema) and/or around (e.g., pleural effusion) the lungs. FIG. 13B illustrates a time 
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domain thoracic impedance signal 1304 and a frequency domain thoracic impedance 
signal 1306 in a subject with pulmonary fluid accumulation in (e.g., pulmonary 
edema) and/or around (e.g., pleural effusion) the lungs. 

As seen in FIGS. 13A - 13B, fluid accumulation can have a profound effect 
5 on the spectral content of a time-domain thoracic impedance signal. Fluid 

accumulation will cause a decrease in the substantially DC component 1302 A and 
1306A (e.g., from C to Z), as illustrated in FIG. 13B) of the frequency domain 
thoracic impedance signal. The fluid accumulation also results in a loss of available 
lung volume for gas exchange. Therefore, patients with fluid accumulation typically 

10 experience rapid, but shallow, breathing. FIGS. 13A - 13B illustrates a shift in the 
higher frequency spectral component 1302B arid 1306B from around 0.17 Hz, for 
the case without fluid accumulation, to around 0.33 Hz, for the case with fluid 
accumulation. The component 1306B also manifests a lower amplitude, for the case 
with fluid accumulation, than the component 1302B, for the case without fluid 

1 5 accumulation. 

FIG. 14 is a schematic diagram, similar to FIG. 2, but illustrating a thoracic 
monitor circuit 112 that obtains a thoracic signal from an acceleration detector 
circuit 1400. The acceleration detector circuit 1400 includes an active or passive 
accelerometer 1402 or like sensor capable of obtaining respiration and/or heart 

20 contraction information from a subject. The accelerometer 1402 includes at least 
one output 1404 that is coupled to at least one input 1406 of an accelerometer 
interface circuit 1408. The accelerometer interface circuit 1408 includes at least one 
output 1410 providing an acceleration-based thoracic signal to the input 206 of the 
analog signal processing circuit 208. One suitable example of the accelerometer 

25 interface circuit 1408, for obtaining an acceleration signal including both respiration 
and heart contraction information, is described in commonly-assigned Kadhiresan et 
al. U.S. Patent No. 5,974,340 entitled APPARATUS AND METHOD FOR 
MONITORING RESPIRATORY FUNCTION IN HEART FAILURE PATIENTS 
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TO DETERMINE EFFICACY OF THERAPY, which is incorporated herein by 
reference in its entirety, including its description of interface circuits for use with an 
accelerometer, such as for obtaining such information. 

It is to be understood that the above description is intended to be illustrative, 
5 and not restrictive. For example, the above-discussed examples may be used in 

combination with each other. Many other embodiments will be apparent to those of 
skill in the art upon reviewing the above description. The scope of the invention 
should, therefore, be determined with reference to the appended claims, along with 
the full scope of equivalents to which such claims are entitled. Moreover, the terms 
10 "first," "second," "third," etc. are used merely as labels, and are not intended to 
impose numeric requirements on their objects. 
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